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Рассматриваются гиперболическое и обратное гауссовское распределения из класса обобщенных гиперболических распределе- 
ний для описания финансовых временных рядов. Предлагается алгоритм оценивания параметров этих распределений с помо- 
щью метода максимального правдоподобия. Апробация алгоритма проведена на примерах эмпирических финансовых времен- 
ных рядов. 


Введение 

Известно, что возвраты большинства финансо- 
вых активов являются лептокуртическими, т. е. 
функция плотности более вытянута в области сред- 
него значения и имеет более тяжелые хвосты, чем у 
нормального распределения [1]. Неудовлетворитель- 
ные результаты прогнозирования, полученные при 
условии нормальности распределения возвратов, за- 
ставляют искать новые распределения и разрабаты- 
вать подходы для обработки эмпирических финан- 
совых данных. Так, МапФеІЪгоІ предложил использо- 
вать устойчивые законы Парето или а-устойчивые 
законы для описания финансовых временных рядов 
[2]. В работах [3, 4] для этих целей было использова- 
но обобщенное /-распределение Стьюдента, в [5] - 
распределение Лапласа. В 1977 г. Вагпсіор- ЫіеЬеп [6] 
описал класс обобщенных гиперболических распре- 
делений (Оепегаіііесі НурегЪоІіс - СН), который стал 
очень популярным в областях теоретической и прак- 
тической статистик. ОН-распределение активно ис- 
пользовалось в физике, биологии и агрономии, а в 
1995 г. ЕЬегІеіп и Кеііег впервые применили его в фи- 
нансах [7]. Указанное распределение имеет ряд 
свойств, которые являются привлекательными для 
описания финансовых временных рядов: 

• СН-распределение позволяет учитывать асимме- 
тричность (известно, что функция плотности воз- 
вратов финансовых активов имеет асимметрию); 


• хвосты СН-распределения тяжелее, чем у нор- 
мального распределения (возникновение ред- 
ких событий, влияющих на форму и вид хво- 
стов, соответствует получению наибольшей 
возможной прибыли или риску наибольшего 
вероятного убытка). 

В статье рассматриваются распределения, кото- 
рые являются подклассами обобщенного гипербо- 
лического распределения: гиперболическое ра- 
спределение ( НурегЬоІіс НУР) и обратное гауссов- 
ское распределение (ТѴоша/ Іпѵеше Оашіап N10). 
Предлагается алгоритм оценивания параметров эт- 
их распределений с использованием метода макси- 
мального правдоподобия. 


Постановка задачи 

Функция плотности обобщенного гиперболи- 
ческого распределения имеет вид: 


ф(х-,Л,р,а,р,8) 
(а 2 - р 2 ) 


(8 2 +(х-/г)) іА_1/2>/2 х 


>/ 2 да х - ІІ2 8 2 К,(8^а 2 -р 2 ) 

хК а-ш (« Ѵ<5 2 +(эг-ц) 2 ) ехр Р(х- /л). 


где \х и 8 - параметры положения и масштаба; 
Р - асимметрии, а - устойчивости. Параметр АеЛ 
характеризует определенный подкласс из семейства 


11 



Известия Томского политехнического университета. 2006. Т. 309. № 6 


00 

А 

Г 01 

Г А-1 

х 

У ехр 


У+- 

0 

ѵ 2 

1 у)) 


обобщенных гиперболических распределений. Для 
хе К функция К х (-) определяется модифицированной 
функцией Бесселя третьего порядка с параметром Я 


сіу, х> 0 


и следующими свойствами: 

1. К, является симметричной относительно Я, т. е. 

К,(х) = К_ х (х). 

1 А 7 -- 

2. Для Я = ±— , имеем К , (х) = , —х 2 е х . 

2 У V 2 

3. К[{х) = --К х {х)-К 1 _,{х). 

X 

4. ЗД = -^ +іѴ 


К х+1 (х) 

К Ах) ' 


5. (Іо ёКАх))' =--К х (х). 
х 

Так, параметр Я=1 приводит к гиперболическо- 
му распределению с плотностью: 

Ъур(х;а,Р,8,/и) = =^х 

28аК і (8 л ]а 2 ~р 2 ) 

хехр(-а^8 2 +(х- /и) 2 + Р(х- /и)). 


Соответственно, при Я = -— получаем обрат- 
ное гауссовское распределение с плотностью: 


К, 


пі ё(х;а,р,8,ц) = 


а8 


ад, 1 + 


х- /и 
8 


71 х/д 2 +(х-/л) 2 

хехр (д-у/а 2 - р 2 + Р(х- /л)). 

Назначение параметров аир для функции 
§Ь(х,Я, у и,а,/3,д) предлагается осуществлять с при- 
менением двух параметризаций [8]: 

я & С^+х 2 

Р = — и а = - . 

8 8 

Тогда функции плотности исследуемых распре- 
делений приобретают следующий вид: 

Ьур(у г А, 8, ц)= , — — х 

2дхД + г ВД) 


Ставится задача оценивания четверки параме- 
тров (г,^,д,/т) гиперболического и обратного гаус- 
совского распределений с использованием метода 
максимального правдоподобия. 

Оценивание параметров 

Предположим, что х„ і=1,...,п - независимые 
наблюдения, тогда функции максимального пра- 
вдоподобия для НУР и N10 -распределений имеют 
следующий вид: 

А = ^\оё(1гур(х 1 ,т,А8, іи)), 

і = 1 

А =^оАпщ(х і ,тА,8,іл)). 

і = 1 

Поиск максимума осуществляется в соответ- 
ствии с выполнением необходимого условия суще- 
ствования экстремума четырех переменных. 

В случае гиперболического распределения: 


(ІЦ _ с 

сіт 8 


сІЦ _ 1 

Ис,~~8 


( и т 

Х(у-а> — г^ х 

ы х/1 + т 2 

хХх/д 2 +(х,.-п) 2 - 
V *•=! 

Г Г ^-А(С)-ѴІ+Ах 


нгд 


С(і+А). 


= 0 ; 


-ид 


+( х і- ѵ) 2 -р) 


= 0 ; 


ац _ с 

дд ^д 2 


сіц _ с 

д/т г 


—ид 

т 


Ч-хЯ+Т" х 


У Я* 1 :' 01 -У>,-<ч 

м Я 5 +(у-р) ;=1 


= 0 ; 


хЯАах 


у-р 


/=1 х/д 2 +(х,. -ц) 2 


= 0 . 


Соответственно, для обратного гауссовского ра- 
спределения: 

г/Я, 2ит тЯ 
— ~ = г~ / х 

дг 1 + г х/1 + т 2 


< 1 ^ 


У - И 


х I А(й,)+^Е(у-р) = °; 

д / о і >1 


ехр 
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X 

+ 


г/Я 2 


1 + ^-|-Ѵі + г 2 


У - Р 


А(у)+7Ё(у -р) = 0; 


. . _ . сѴі+А 

пі8(х,г,С,д,п) = х 
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Естественные науки 
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возвраты, т. е. 1п — — , где 5 Л - значение котировки 
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валютной пары в день п. 

В табл. 1 сведены основные статистические ха- 
рактеристики по всем валютным парам, в табл. 2 
показаны оценки параметров для НУР- и N10 -ра- 
спределений. 


Алгоритмы решений двух указанных систем не- 
линейных уравнений относительно неизвестных 
параметров (т,^,8,рі) в предположении единствен- 
ности экстремума осуществляются на базе рекур- 
рентных процедур методом наискорейшего спуска, 
сопряженных градиентов и т. п. 

Для N10 -распределения в качестве начальных 
значений искомых параметров используются ре- 
шения, определяемые следующей системой ура- 
внений: 


Е(Х) = /и + 8т; 


Ѵ(Х) 


5 2 (1+т 2 ) 

’ 


8{Х) = 3 


К(Х) = - 

С 


л/с ѵг 


1 + 4- 


+ г~ 

Л 


1 + г 


Здесь Е(Х) - математическое ожидание выбор- 
ки, Ѵ(Х) - дисперсия, 8(Х) - коэффициент асим- 
метрии, К(Х) - куртозис. Решение этой системы 
уравнений приведет к нахождению параметров 
(г, і^,8, /и), которые могут быть использованы в ка- 
честве начальных значений для оценивания пара- 
метров гиперболического и обратного гауссовского 
распределений. 

Соответственно, для НУР-распределения в ка- 
честве начальных значений используются данные 
НурегЬоІісБіМ К-раска§е$ [9]. 


Таблица 1. Основные статистические характеристики 


Котировки валют 

М{Х) 

Ѵ(Х) 

5(Х) 

К(Х) 

ЕІІК/1150 

-0,0001 

0,000029 

0,21673 

0,09510 

СВР/ІІБО 

-0,00002 

0,000025 

0,23685 

0,03847 

ІІБЦ/ІРѴ 

0,00031 

0,000028 

-0,43221 

1,03296 

ІІБЦ/СНЕ 

0,00019 

0,000035 

-0,36500 

0,22951 


Таблица 2. Оценки параметров для НУР и N10 -распределений 


Котировки 

валют 

Параметры 

т 1 С 1 <5 I Р 

і НѴР-распределение | 

ЕІІК/ІІБЦ 

0,34582 

20,54527 

0,02205 

-0,00832 

СВР/ІІБЦ 

3,66805 

16,96968 

0,12459 

-0,06550 

иБй/іРУ 

-0,28113 

5,81484 

0,01085 

0,00417 

иБР/СНЕ 

-0,43780 

9,93258 

0,01571 

0,00813 

I МІС-распределение | 

еіік/іібц 

0,07070 

7,07653 

0,01420 

-0,00114 

СВР/УБЦ 

5,25364 

10,64734 

0,12843 

-0,04730 

иБй/іРУ 

-0,20705 

4,60074 

0,01119 

0,00277 

1150/СНЕ 

-0,28605 

2,11668 

0,01644 

0,01459 


Из табл. 1 видно, что распределение логариф- 
мических возвратов эмпирических данных харак- 
теризуется асимметрией, наличием куртозиса и не 
может быть описано распределением Гаусса. Ис- 
пользуя оценки параметров из табл. 2, можно по- 
строить гиперболическое или обратное гауссовское 
распределения, позволяющие наиболее адекватно 
описать эмпирические финансовые данные. 

Выводы 


Эконометрических анализ данных 

В работе используются данные по валютным 
парам ЕШ/ШБ, СВР/ШБ, ШБДРУ, ШБ/СНЕ 
за период с 03.01.2004 по 29.09.2006 г. (всего 455 
значений) в качестве эмпирических временных ря- 
дов для оценивания параметров НУР- и №С-ра- 
спределений. Данные состоят из дневных цен зак- 
рытия, из которых формируются логарифмические 


Рассмотрены гиперболическое и обратное гаус- 
совское распределения из класса обобщенных ги- 
перболических распределений и предложен алго- 
ритм оценивания параметров этих распределений с 
помощью метода максимального правдоподобия. 
Полученные оценки параметров используются для 
описания финансовых временных рядов и постро- 
ения адекватных математических моделей. Алгоритм 
успешно апробирован на эмпирических данных. 
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